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Abstract 

The calculation of the finite lattice spacing corrections for 1=2 vrvr scattering is carried out 
for isotropic and anisotropic Wilson lattice actions. Pion masses and decay constants are also 
determined in this context. These results correct the phase shift calculated from the lattice, which 
is connected to the scattering length and effective range in this low energy scattering process. 
When in terms of the lattice-physical parameters for either Wilson action, these lattice spacing 
effects first appear at the next-to-leading order counter-terms. 
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I. INTRODUCTION 



Numerical scattering calculations in lattice QCD are being performed by several collab- 
orations. These calculations are performed through the analysis of two hadrons in finite 
volume [l|, Q, H, S] • One such scattering that has gained much attention in the field is 1=2 
TTTT scattering. Such numerical calculations (usually involvin g p hase shifts and scatte ring 
len gth s] have been calculated using Wilson lattice actions [5|, |6|, 

UMM B El, [H, E, H 



along with a several other lattice actions [21|, |22|, |23|. Additionally, 
there are currently only two full y d ynamical, 2+1 flavor calculations of 1=2 vrvr scattering, 
which use mixed lattice actions I 



Regardless of the action, unphysical lattice arti- 
facts due to the finite lattice spacings exist in the numerical results of these calculations. 
Therefore, measures should be taken in order to remove these effects so that the results from 
the lattice can best represent the continuum limit. The analysis in this paper is applicable 
for both isotropic and anisotropic Wilson actions. 

Effective field theory (EFT) provides a framework by which one can remove these un- 
physical effects. Lattice spacing effects were first made explicit in chiral perturbation theory 
(xPT ) by Sharpe and Singleton [26| . For the Wilson action, the chiral breaking terms that 
depend on the lattice spacing can be accounted for in a similar way to the chiral breaking 
quark mass. Such methods have been extended to mixed-action, partially quenched theories 
for mesons through 0{a) and 0(^) ^ 



28, 29 



and 0{a^) 

\Mm, S3. Additionally, Ref. IM carries out multiple meson scattering calculations (in- 



^ and baryons through 0{a) \3i 

13 . Additionally, Ref. ^ 
eluding 1=2 tttt scattering) for mixed lattice actions and shows that for actions with chiral 



valence fermions, mesonic scattering parameters in terms of the lattice-physical parameters 
will have no counter-terms dependent on lattice spacing through next-to-leading order. Al- 
ternatively, this work calculates these lattice spacing effects^ for 1=2 vrvr scattering for the 
chiral breaking Wilson fermions in both valance and sea sectors, and shows that these finite 
lattice spacing effects first appear in the next-to-leading order counter-terms for this action. 

Many collaborations are now using anisotropic lattices (lattices with different temporal 
and spacial lattice spacings) as opposed the the usual isotropic lattices. Such lattices can 
probe higher energy states (inverse time spacings a^^ ~ 6 GeV) and allow for a greater 
resolution (more data points). However, anisotropic lattices lead to new lattice artifacts, 
including terms that explicitly break hypercubic symmetry. Recent work has derived these 
anisotropic lattice artifacts in xPT for 0{a) and O(a^) for mesons and baryons 35|]. There 
are several 1=2 vrvr scattering results published for anisotropic Wilson lattices [lo|, 17, 20|, 
which can benefit from removing these additional lattice artifacts. 

This paper presents the 1=2 tttt scattering results from the isotropic xPT and the 
anisotropic xPT . The pion mass and decay constant are also determined in this context. 
Sec. ITT] presents scattering on the lattice and defines the relevant quantity, k cot So, used to 
make comparisons between xPT and the actual lattice calculation. Next, in Sec. IIII A[ the 
continuum scattering theory from xPT is formulated in the context of this paper (originally 
worked out before 36, 37, H, 39|). Then, in Sec. IIIIBl the continuum result is extended for 
the isotropic Wilson lattice and finally, in Sec. lIII CI the result is extended for the anisotropic 
Wilson lattice. 



^ These effects are in addition to the finite volume corrections to 1=2 tttt scattering from Ref. [34 
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II. SCATTERING ON THE LATTICE 



The Euclidean two-hadron correlation function in infinite volume gives no information 
about the Minkowski scattering amplitude (except at kinematic thresholds) j4o|. However, 
when the correlation functions of two hadrons in a finite box are analyzed, the resulting 
energy levels are given by the sum of the energies of these two hadrons plus an additional 
energy of interaction, Ai?, which is related to the scattering phase shift, 5/ [H, 0, H, 0] ■ The 
/ subscript here represents the partial wave contribution of the phase shift. In the infinite 
volume, the relation between the total scattering amplitude, T{s,9), and the partial waves 
amplitude, is given by 

oo 

T(s,^) = ^(2Z + l)P,(cos^)tKs), (1) 

1=0 

where s = 4(m^^ + k"^), and k is the magnitude of the 3-momentum of the incoming particle 
in the center-of-mass frame. The partial scattering amplitude t;(s) is related to the phase 
shift, 6 1 by 



ti{s) = 3271 J ^ --[e^^-^'^^) - 1] = 32nJ ^ ^ . (2) 

^' \l s-Am^^2i^ ^ y s-im^^ cot 6i-i 

These relations allow one to compare the calculated scattering amplitude (in xPT ) 
to the lattice calculation of 6i. The s-wave (/ = 0) scattering amplitude is the dominant 
contribution to the total scattering amplitude in most low energy scattering processes and 
gives the cleanest signal in the lattice calculation. The s-wave projection of the continuum 
scattering amplitude, to{s), is 

1 



to{s) = -J T{s,e) dicosO). (3) 

This s-wave scattering amplitude will be the scattering amplitude analyzed throughout 
the rest of this paper. Following the discussion in Ref. through one loop order in 



perturbation theory in Minkowski space, to("S) can be written as 



T^O{S) -Iq {S) +Iq (S) ^IIq [S) - _ ANLO,R), . _ -ANLO,!), x' 



where t^Q^\s) is the leading order s-wave scattering amplutide, and t\^^^'^\s) {tl^^^'^\s)) 
is the real (imaginary) part of the next-to-leading order s-wave scattering amplitude. At 
this point, it is advantageous to introduce a i^-matrix, which is defined through one loop as 

K(s) = («°'°'(^')' (5) 

Taking the real part of the reciprocal of Eq. ([2]) and Eq. (jlj) and comparing to Eq. (Q, 
one gets the relation 



K{s) \to{s) J 327r 



HT7^^] ^cot6o{s), (6) 
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where 



RelJ-]^ ^^—d-'l^^] (7) 

W-l)J (flei„(s))'+(/mi„(s))" 4^°'MV 4^°'W / 

It is worth noting that when keeping terms though one loop, Im to(s) does not contribute 
(it contributes at the next order). Combining Eq. ([5]), Eq. ([6]), and Eq. ([7]), one arrives at 
the continuum resuh 

/ 1 \ 1 / t^^^^'^^(s)\ 

*cot.„(.) ^ 16.v/;«e(^) . - ^^). (8) 

As previously mentioned, lattice scattering calculations are performed in Euclidean space 
at finite volume. The Euclidean amputated four-point correlator from the lattice, ro(s), is 
given by 

where Ato(s) represents all of the non-physical lattice artifacts (lattice spacing errors, finite 
volume errors ,etc.), s is related to the energy of interaction, AE, and iS is a universal 
function of s jsl, |4l|, |42| . If both pions in the box start with no external momentum, then 



s = {AE + 2m^f. In this paper, the only effect from lattice artifacts that will be included 
in Ato is the lattice spacing effect. Manipulating Eq. ([9]): 



^ ' ' Afo(^) 1 ^/^ (s-4m^2)L2 



167r-\/i 

" fccot5o(3) - - ^^((^^^ 

The energy states are given by the poles of Eq. ffTOl) . which are given by [3. 

ttL V 167r 



(10) 



fccot5o + A(fccot5o) = — 5( ^ T7r^^)> (11) 



where 

A(fccot(5o) = -167rv^ ^/°y . (12) 

In general, most lattice calculations give their results in terms of the scattering length, 
a^^^. One can extract the scattering length and the effective range, t^^^, via the expansion 
of fccot 5q: 

kcoi5o = ^ + \riZ^k^ + --- . (13) 

It is important to note that the prescription given above for finding the scattering length 
and effective range implies that the lattice artifacts, A(fccot5o); have already been sub- 
tracted before the expansion. In this paper, continuum results are given in terms of k cot 5q 
and lattice artifacts are given in terms of A(A;cot(5o). 
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For results given in terms of the scattering length and effective range, one can relate 
Eq. (fT3|) to the left hand side of Eq. (fTTj) to arrive at 



k cot 5o + A(A: cot 6,) = {J^^ + A (J^) ^ + \ [rif + Ar^f ) fc^ + ■ ■ ■ , (14) 

where 

a(-^) =A(A; cot 5o) 1^2=0, (15) 



and 



_ o ^(A(fecot(5o)) 



(16) 

While these relations are not too complicated, they do add additional steps to the calcu- 
lation when compared to working with only k cot 5o and A[k cot 5o). Therefore, if one wants 
to extract the scattering length and effective range from the lattice calculation, one should 
first subtract A(fc cot 5q) from the right hand side of Eq. (|TT1) and then expand to determine 
the individual parameters^. This paper relates fccot(5o and A(A;cot5o) to the effective field 
theory of the lattice. 



III. CHIRAL PERTURBATION THEORY RESULTS FOR fccot^o AND A(/i:cot(5o) 

In leading order (LO) and next-to-leading-order (NLO) chiral perturbation theory 
(xPT ), it is necessary to introduce several undetermined low energy constants (LECs) 
in order properly account for corrections and counter-terms. The number of independent 
LECs in the continuum depends on whether there are two flavors or more. For 1=2 tttt 
scattering being calculated here, only two flavor xPT {SU{2)l ® SU{2)r chiral symmetry) 
is needed for extrapolation. 



A. Continuum 



From the continuum xPT Lagrangian, one can predict numerous results for different low 
energy processes involving hadrons. However, the extent of the accuracy of these predictions 
are ultimately tied to how well the LECs are known. For this reason, there has been much 
effort in the lattice community to try to determine these values (isl . 44. 45. 



46, 47 



The continuum Lagrangian in xPT is determined order by order in Bruq and p . The 
Lagrangian through 0{p'^) for two flavors is given by [sT" 



^ Current numerical calculations can only determine k cot So for a limited number of k values. This leads 
to inaccuracies in expansions of fc^ and adds difficulty to finding the effective range. 



5 



f2 TD £2 



(17) 



where / ~ 132 MeV, are the original Gasser-Leutwyler coefficients defined in Ref. [3 
and 



E = exp 



2i0 

T 



_7ro_ 

7r 



TT 



V2 



m 
m 



At LO, the resulting condensate is 



B = lim 



(18) 



(19) 



From Eq. ( fT7|) . one can calculate the physical values for the mass of the pion (m^r) and 
the pion decay constant {f-^) to LO and NLO. These expressions are given by 



1 tHj 



U = f 



2 n 



2 777; 

1 - --^il{m^) + 24—^ 



(20) 



(21) 



where and I{mT^) are defined below in Eq. fl23|) . When evaluating the scattering am- 
plitude from xPT , one has the option of either expressing the answer in terms of the bare 
parameters (/ and m) or in terms of lattice-physical parameters (/j^ and m.,^). Throughout 
this paper, the bare parameters will always be eliminated from the scatting amplitude. The 
continuum 1=2 tttt scatting length at infinite volume is given by 



T, 



cont 



2 



2(3s - 4m^2 



(s — 2m7r 



2\2 



2)S - 2m^^ ' \ r 2 " ' i'^i'T^Tv) H 71 iJirn-^.ps 



3// 



+ 



2 4\ 



t{t — s) — 2tmT,^ + AsniT^'^ — 2'mT,'^~\iJ'{m^,pt) 



3// 

-|- zTT^ [3(-u^ — m^*^) -|- u{u — s) — 2um.„'^ + ism.,,'^ — 2m^^]ij7'(m^,pu) 



3// 



1 



■ [2s^ + 6 



9(471^)2 

^[(t-2m^y + {u-2m^Y] 

JTT 

2i2 



[2{s - 2m^^)^ + {t- 2m^2)2 + {u - 2m^'^f] 



-■i-pr + 71 

/tt /tt 



(22) 
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where 



m 



2Bm 



d^k 1 



(23) 



This scatting amphtude includes all the partial wave contributions (this T is the same 
as the T{s, 9) that appears in Eq. ([1])). When projecting on the s-wave, expanding through 
0{k'^ /rriT^'^) and using Eq. (EL the result for A;cot5o(for the regularization and renormaliza- 
tion scheme defined in Ref. [38j) is 



k cot 5q 



1 - 



(4vrA 



3 In 



^ (3 1 

2 V (47r/.)2 



11 m 

L 3 



+ ■ 



(24) 



where ^"^(/i) and ^^^(/i) are linear combinations of the Gasser-Leutwyler coefficients given 
by H 



''TTTT 



(^) = - 4(47r)2 (4(£f (/x) + £f (^)) + (£f (^) - £f (/i)) 
(/i) =4(47r)2(l2£f (/i) + 4£^(^) + 7£f (/i) - 3£f (/i)). 



(25) 



The superscript R represents the renormalized Gasser-Leutwyler coefficients with scale- 
dependence. 

In order for these continuum predictions from xPT to be useful in a physical context, 
one needs to determine the values for ^"^(/i) and £^^{fi), which are undetermined from 
yPT alone. While numerous values have been quoted for the Gasser-Leutwyler coefficients 
[47|, it is still beneficial to determine these values with more precesion. Therefore, it is 
prudent to use various lattice calculations at different pion masses determine these values. 
However, since the lattice calculations are performed with discretized space and time, one 
needs to remove these lattice artifacts to extract the continuum result. 



B. Isotropic Discretization 

To calculate finite lattice spacing corrections to 1=2 vrvr scattering for the isotropic Wilson 
action, one can follow the same steps done in the continuum case, but starting from a 
Lagrangian which includes these lattice spacing artifacts. The analysis on lattice spacing 
effects was done for the Symanzik action by Sheikholeslami and Wohlert J48[| . From this 
analysis, the Lagrangian was made explicit in xPT by Sharpe and Singleton |26[], followed by 
Bar, Rupak, and Shoresh [23, |28l]. The power counting they used for this xPT Lagrangian 
is 

asW ~ Bmg ~ ~ e, (26) 
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where as is the lattice spacing (same spacing in space and time direction) and is a 
condensate defined below. The simplified two fiavor Lagrangian to 0{e^) is 



^iso l~'Cont ~l~ ^(~'iso- (^'^) 

A£,„ =^^tr(Et + E) + <'"^ + '°f-'"^^° tr(.n,Et + E™,)tr(Et + E) 

+ ^ j,,] 2 ^ !^:i^tr(9„ES^Et)tr(Et + E) (28) 

This Lagrangian is similar to Eq. (IT7|) with one new term at LO and three new terms at 
NLO. All new terms are proportional to or a^, which will vanish in the continuum limit 
as as ^ 0. At LO, there is a new condensate given by 

W = hm csw • (29) 

The new LECs at NLO are given by w^, w'^, and W4. All these new terms obey the same 
symmetries as the Lagrangian in the continuum case and break chiral symmetry in a similar 
way to the quark mass. It should also be noted that these new LECs depend on as In as 
well (as opposed to the Gasser-Leutwyler coefficients that have no dependence on the mass 
term) . 

Furthermore, Ref. 4^ showed that the axial current (needed to calculate /^r) has an 



additional term at this order in xPT given by 

AA;: = 2awA9^tr(a'^(S-St)). (30) 

This term (which can also be inferred from Ref. [50|) leads to modifications of the LECs 
as well as the coefficient in front of the chiral logarithm in f„. Thus, this coefficient has 
dependence on the lattice artifacts at NNLO. Ref. [49| also points out that the condition 
for fixing the renormalization factor, Za, of the lattice currents needs to be mapped onto 
xPT ^. From this Lagrangian, and m-Tr through NLO are 



_ ^^2 ^ 2asW) + -^[Am^HX{m^) - {vr? + 1asW)%Z{m^)\ 

Jtt J TV 



U = f 



2 , , ^ asWi 

Jtt Jtt Jtt 



(31) 



(32) 



where Weff in /^r includes W4 and wa and can vary based on the given renormalization 
condition for the axial current. 



The condition for fixing Za is chosen by individual lattice calculations. Ref. [49| shows an example of this 
in xPT by using the chiral Ward identities in infinite volume, which leads to /jr being free lattice artifacts 
until NNLO. 
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To acquire the continuum result of these quantities, one needs to remove all the terms 
with dependence on or a^. The resulting 1=2 vrvr scatting amplitude with the physical 
parameters restored is given by 

where the rriT, and in Tcont are given by Eq. (l3Ti) and Eq. (l32i) and ATj^o is given by 



ATiso = --f^) - 16(1^3 - 24) — 32(^3 -W3 + £3) 

Jtt L Jit 



+ 8K// - ^4) ^ |- (34) 

It is worth noting that by restoring the physical parameters, does not appear, and 
QgW only appears in the terms containing the LECs (the i and w terms). This is a bit 
different than the continuum case where one could eliminate with only m^^. Now, one 
eliminates with (m^^ — 2asW), and thus, several LECs are multiplied by factors of asW. 
In addition, all of the continuum results without LECs remain unchanged since each vertex 
will only contribute m^r^ when the physical parameters are restored. Using the relation in 
Eq. ffT^ . the resulting artifact for the isotropic Wilson lattice, A(/c cot 5o)jso is given by 



A A;cot5oU~77^< (<M — 
27r V "^TT 



.^^M^ + 7<M- 4 / 2 
2 V J m ^ 



- <^(/x)-^ + 7<';(/x)^^-^ — ^ + --- \, (35) 



where 



KM = - m^rKif^) - <fM) - (/.) + (/.)), 

= - 16(47r)2(«;f (^) - wfi^^) + £f (/i)), 
(/.) = - 8(47r)2 (7<(/x) - 3<//i) - l<(/i) + 3£f (/.)) . (36) 



As seen in the results, the artifacts from the isotropic lattice that are present in the final 
form are either linear or quadratic in a^. By using results that differ in lattice spacing, 
one can pick off the coefficients of these artifacts and remove them from the final result. If 
one is working with a perfectly clover-improved Wilson lattice, this would remove all 0{as) 
effects leaving only the 0{a1) effects. It is also important to note that there is no physical 
information gained by determining specific LECs that are a result of the isotropic lattice 
spacing (the individual w terms) unlike determining specific Gasser-Leutwyler coefficents. 
Therefore, the useful coefficient to extract is the linear combination of these terms so they 
can be removed from the final result. 

At this point, one can compare these lattice spacing effects for the Wilson action to 



those found for the mixed action (with chiral valence fermions) in Ref. 29, 51 . For this 



mixed action case, when in terms of the lattice-physical parameters, there is no lattice 
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spacing dependence through the NLO counter-terms. In contrast, when the lattice-physical 
parameters are restored in the Wilson action, these effects first appear at the NLO counter- 
terms. Thus, these additional effects that are not present at NLO in the mixed action 
calculation will need to be removed for the Wilson action calculation. 



C. Anisotropic Discretization 

With the recent formulation of xPT for the anisotropic lattice [s^ , one can begin calcu- 
lating corrections to various quantities of interest on the lattice. This process is, in general, 
very similar to calculations in the continuum and isotropic lattices, but one now picks up 
additional terms that are a result of having different spacial and temporal spacings. To help 
extract these effects in a more simplistic notation, the superscript ^ has been added to all 
the new terms resulting from this anisotropy. In practice, the anisotropic lattice picks up 
two new non-perturbative parameters: the parameter ^ = as/ at which is the measure of 



anisotropy and the parameter z/, which is used to correct the "speed of light" [52|, |53|, [54 
By setting both parameters to 1, the isotropic limit is recovered. In addition to the W 
condensate defined in Eq. fl2^ . we pick up a condensate that is given by 

W^= lim4^(^^)^(^^). ^''^^^„r'^ - (37) 

where is a vector that breaks hypercubic invariance. It is important to note that the con- 
vention chosen for appears in the anisotropic xPT and its observables. For convenience, 
we choose this vector to be = (1,0). The condensates and the anisotropic paramaters 
are related at the classic level by (with Wilson coefficients Ts = = 1) 

W oc Csw oc z^, (38) 



2 

In the isotropic limit when ^ and v are set to 1, the isotropic condensate will remain and 
the anisotropic condensate will vanish. Using a similar notation throughout, all terms that 
appear with a ^ superscript will vanish when as = at and u = 1. 

The power counting convention used in Eq. (12^ for the anisotropic Lagrangian is 

a.W ~ a, ^ e. (40) 

Writing this new Lagrangian in the form of Eq. fl27p . the two-flavor anisotropic 
xPT Lagrangian through 0{e^) is 

^aniso ^cont ~l~ ^^iso ~l~ ^^aniso- (^-^) 

where 



+ 4 



^-f— tr(9^Sa'^St)tr(St + S) + ^ f M^M,(.9^Sg'^St)tr(St + S). (42) 
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In addition to the anisotropic condensate mentioned above at LO, there are five new 
LECs at NLO as a result of this anisotropy. Four of these new LECs obey the same sym- 
metry structure as the isotropic terms, however the a;f term additionally breaks hypercubic 
invariance. Therefore, this term only corrects the time derivative, but not the spacial one 
(for the convention of chosen). As a result, fj^ is parameterized by two constants; 
which is measured in time, and which is f-„ measured in space. This leads to one cor- 
rection for the space-measured and a separate correction for the The pion mass (771^^), 
the time-measured pion decay constant (/*), and the space-measured pion decay constant 
(/^) through NLO are 



[m 



+ 2asW + 2asW^) + — [Am^'^ilim^) - {ir? + 2a,W + 2a,W^)iX{m^)\ 



3/ 



a^Wm 



^W^— 7^ h IQWn TT- 

J n J IT 



tUsKK -III, ^i:(asW^Y 



J TV 



J TT 



J TT 



J n 

fS 



r ^ Tft^ Qj 



ft 



m 



ft- 



ft- 



f 1 - j^^A^tt) + + 4«;4^ + M 



a 



ft- 

J TT 



ft- 

J TT 



J TT 



(43) 
(44) 
(45) 



where, as in the isotropic case, the u^e// depends on the renormalization condition for the 
axial current (this LEG is the same for both /* and /^). 

For the rest of this section, all calculations are given in terms of m^r and As mentioned 
before, how one accounts for the effect of this hypercubic breaking term depends on the 
convention. In the convention used here, only /* sees the effect of this term and does 
not. The scatting amplitude is given by 



where the m.,, and /* in Tcont are given by Eq. fj43l) and Eq. 



(46) 

and /STaniso is given by 



^Taniso = "(y^l " 16(^1 - 2h)—^2 32(w;| - w| + £3) 



2 



- 32{wl -W3-wl + 2i, 



{asW){a,W^) 



3 



Most of the effects in this scattering amplitude are similar to the isotropic case, except 
now there are also expansions in agW^ in addition to the expansions in asW. Thus, as 
expected, if all the anisotropic effects are removed, only the isotropic limit remains. The 
only new symmetry breaking effect is the Wi term which is not a hypercubic invariant term. 
In other words, all of the hypercubic breaking due to anisotropy at this order is contained in 
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this term. However, it's effects in ATaniso appear as just another contribution to the hnear 
combination of the LECs in front of the term asW^. Therefore, it is difficuh to determine 
the effect of the the hypercubic breaking term alone from ATaniso since its effects will be 
mixed in with the other anisotropic LECs^. The resulting artifacts for the anisotropic Wilson 
lattice are the isotropic artifacts, A{k cot 5o)iso) and the anisotropic artifacts, A{k cot 5o)amso- 
Therefore, the total effect of the lattice artifacts due to lattice spacings are 

A{k cot So) = A{k cot 5o)iso + A(A; cot 5o)amso- (48) 
The anisotropic lattice artifacts are given by 



(49) 



where 



^i^'M = - 8(47r)2(«;f (/.) - wif^{f,) + (/.) - 2ii{i,) + £«(/.)), 

= - 16(4vr)2(i.f (/.) - «;3^(/x) + £3^(/.)), 
wi^M = - 16(47r)2(^f (^) - wfi^,) - (/.) + 2£f (/.)), 

^i;(/^) = - 8(47r)2 (7^f (/.) + 9wiff{i,) + ^f^(/.) - 14£f (/.) - <(/.)) . (50) 

As a result, the ultimate effects of the anisotropic lattice on A(A;cot5o) are more terms 
that will require variation of and at independently to fit. Three more terms require fitting 
to correct the constant term in the expansion and one more term requires fitting for the 
0{k'^ /mT^) term in order to remove it's effects. As mentioned with the isotropic correction, 
no physical information is gained by picking off the individual anisotropic LECs. Analogous 
to the isotropic case, the anisotropic Wilson action first has lattice spacing dependence at 
the NLO counter-terms. The aim here is to determine these linear combinations and remove 
their effects from the lattice measurements. 



IV. DISCUSSION 

In this work, 1=2 tttt scattering was calculated from the xPT for isotropic and anisotropic 
lattice spacings. Also, connections between these xPT calculations and the k cot value 
measured from lattice calculations were illustrated. When A(A: cot 5q) is given in terms of the 
lattice-physical parameters, these lattice spacing effects first appear at the NLO LECs and 
can be removed from the result of the lattice calculation. However, A{k cot 5q) has numerous 
undetermined linear combinations of LECs, which would need to be determined (by fitting 



The total hypercubic breaking effect would be more visible from the differences of /* and 



12 



several different lattice spacings) in order to successfully remove it from the lattice result. 
Therefore, as more lattice calculations of vrvr scattering are completed (for both isotropic and 
anisotropic lattice spacings), these combinations of LECs can be determined better, which 
will result in a more accurate result after these lattice artifacts are removed. 
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